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UNIT III: HIGHER-ORDER LINEAR EQUATIONS 



We give a comprehensive development of the theory of linear differential equations with con- 
stant coefficients. We use the operator calculus to deduce the existence and uniqueness. We 
presents techniques for finding a complete solution of the inhomogeneous equation from solu- 
tions of the homogeneous equation. We also give qualitative results on asymptotic stability. 



f(t) is a continuous function on an interval I. The letter L stands for the (homogeneous) differen- 
tial operator. It is easy to see that L : C n (I) — > C(I) is linear, where C k (I) is the space of functions 
differentiable k times on /. 

As for the second-order equations treated in Unit II, the principle of superposition and the principle 
of the complementary solution apply to (11.1). 

Principle of Superposition. If Lu = and Lv = 0, where L is given in (11.1), then L(c±u + C2V ) = 
for any constants c\ and C2- 

Principle of the Complementary Solution. If u is a particular solution of Lu = f, where L is 
given in (11.1), and if v is any solution of Lv = 0, then L{u + v ) = / and every solution of Ly = f 
can be obtained this way. 

Therefore, the general solution of (11.1) is given as 

y = y P + Vh, 

where y p is a particular solution of (11.1) and yt is a solution of the corresponding homogeneous 
equation 



The characteristic polynomial. We try y(t) = e xt , X G C, as a solution of the homogeneous equa- 
tion (11.2). Since (e xt ) = \ k e xt , the substitution yields 



Moreover, since e xt is never zero, Le = if and only if A is a root of the characteristic polynomial 
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(11.2) 



Ly = y^' + Pl y {n l > + • • • + p n _ iy ' + p n y = 0. 



Le xt = (A™ +p\\ + • • • +p n . 1 \ + Pn )e xt = 0. 



(11.3) 



PL (X) = X n +p 1 X n - 1 + 



+ p n -i\+p n . 



of L. 



l 



Example 11.1. We recall the results for the second-order equation 

(11.4) Ly = y"+py' + qy, 

where p , q are constants. The roots of the characteristic polynomial are 

- P ± Va 2 

A= , A=p -4q. 

If p 2 > 4q then y = e xt with the above A are solutions of Ly = 0. 

We will show how the exponential substitution y = e xt applies to solve the general differential 
equations of all order (11.2). 

The operator calculus. The study of linear differential equations become easier if we introduce 
an abstract symbol D = d/dt for the operation of differentiation. A word of caution. The symbol 
d is used for differentials, e.g., D(t 3 ) = 3t 2 , but d(t 3 ) = 3t 2 dt. 
As an operator, D is linear. That is, 

D(u + v) = Du + Dv, D{cu) = cDu 

for all differentiable functions u, v and for any constant c. 
We now list some properties of D. By definition, 

D° = id, D k = ^-r A: = 1,2,.... 
dt k 

Moreover, 

(11.5) DW k = D j+k , (D j ) k = D jk , j, k = 1, 2, ... . 

The proof is left as an exercise. 

With the notation D, we may write the differential operator in (11.2) as 

Ly = (D n + p 1 D"- 1 + • • • + pn^D + p n )y. 

Then, it can be recognized that the first factor of the right side is pl(D), the characteristic polyno- 
mial pl evaluated, formally, at D. In this sense, we say L = pi{D). 

Linear operators with constant coefficients are permutable, in the sense of the following. 

Lemma 11.2. Ifp(D) = o,jD^ and q(D) = ^ bkD k are two linear differential operator, where cij , bk 
are constants, then, 

p(D)q{D) = q{D)p{D) = ^^jh^ +k . 
The proof uses (11.5) and it is left as an exercise. 

Remark 11.3. The above lemma is not true of linear operators with variable coefficients. For exam- 
ple, 

D(tf) = (if)' = tf + f = (tD + id)f, 
where / is a differentiable function of t. In other words, Dt = tD + id. 

In many applications, one takes trial solutions of the form e x tu where A G C and u is a function 
with a certain degree of smoothness. Thus, it is useful to know how such a function works with 
the operator D. 

Lemma 11.4 (The Exponential Shift Laws). If p is a polynomial and A is a constant, then 

p{D){e xt f) = e xt p{D + \)f. 
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Proof. By the rule for differentiation, 

D{e xt f) = e xt Df + Xe xt = e xt {D + X)f. 

Hence, 

(D-a)(e xt f) = e xt (D-a + X)f 
for any constant a. Then, by induction, 

(11.6) (D - a) k (e xt f) = e xt (D -a + \) k f. 

Finally, by the fundamental theorem of calculus, any polynomial p can be factorized as 

p(D) = (D — ai ) k '(D - a 2 ) k2 .-.(£>- a m ) fc "% 

where aj G C are roots of the polynomial and kj ^ 1 are the corresponding multiplicity. The 
assertion then follows by Lemma 11.2 and (11.6). □ 

As a consequence, moreover, we have 

(D - X)(e xt f) = e xt Df, (D - \) k {e xt f) = e xt D k f. 

Exercise. If a is not a root of the polynomial p , then show that 

b{t) 



e at 



p(a) 

is a particular solution of the differential equation p{D)y = e at . 

We now present our main result. 

Theorem 11.5. If A* is a (complex) root of multiplicity k of the characteristic polynomial p(X) = X n + 
PiX" 1 ^ 1 + • • • + Pn-i^ + Pn of the linear differential operator p(D) with constant coefficients, then the 
functions i r e A **, where r = 0, 1, . . . , k — 1, are solution ofp(D)y = 0. 

Proof. By the exponential shift law, it follows that 

(D - A*) fe (fe A **) = e x " t D k f = 

for r = 0, 1, . . . , k — 1. 

On the other hand, p(X) must contain the factor (A — A*) fc , and hence 

p(D) = (D — X*) k q(D), q{D) = \{ (D - X,p . 

Finally, by Lemma 11.2 

p{D)(fe Kt ) = q{D){D - A*) ^V**) = 0. 
This completes the proof. □ 
Corollary 11.6. If 

p(A) = (A-A 1 ) fcl (A-A 2 ) fc2 -(A-A m ) fc ™ 

then the functions t T e x ^ , where r = 0, 1, 2, . . . , Aj_i and j = 1, 2, . . . ,m, are solutions of the differential 
equation p(D)y = 0. 
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Complex solutions. An interesting feature of the analysis via the operator calculus is that many 
problems are best solved by the use of complex-valued functions, even when the coefficients are 
real, e.g. (11.4) whenp 2 < 4q. 

Theorem 11.5 holds whether the coefficients pj of p{D)y = are real or complex. In fact, al- 
though t is interpreted as real (in particular in the discussion of stability) the operator calculus 
and the solutions constructed with it apply equally well to functions of the complex variable. But, 
when the coefficients pj are real, then we may obtain a stronger conclusion. 

We recall that the complex roots of a polynomial come in pair, fj, ± iu, n,u e K. We also recall 
the complex exponential e^ ±lv = e^(cos u ± iu). 

Lemma 11.7 (Principle of Equating Real Parts). If a complex-valued function y(t) = u(t) + iv(t), where 
u and v are real-valued functions, satisfies the differential equation (11.2) with real coefficients, then u(t) 
and v(t), the real and imaginary parts of y, both satisfy (11.2). 

Proof. Let Ly = Liu + iv) = 0. Since the coefficients of L are real, taking the complex conjugate 
we have Ly = L(u — iv) = 0. Then, by linearity, 



It holds true for differential equations with real-variable coefficients. 

Corollary 11.8. Each pair of complex roots fi±iv, fj,,v G M, of multiplicity k of the polynomial p gives 
real solutions t r e ut cos vt , t r e ut sin vt, where r = 0, 1, . . . , k — 1, of the differential equation p(D)y = 0. 



u = 



y + y 

2 



and v 



y-y 

2i 



both satisfy (11.2). 



□ 
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